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NONLINEAR  STABILIZATION  OF  THE  FARLEY-BUNEMAN 
INSTABILITY  BY  STRONG  E  X  B  TURBULENCE 

It  is  well  known  that  in  the  absence  of  a  magnetic  field  the  two  stream 
instability  can  occur  in  a  homogeneous  plasma  when  the  electron  drift  velo¬ 
city  with  respect  to  the  ions  exceeds  the  electron  thermal  velocity.1 
Farley2  and  Buneman2  have  shown  that,  in  the  presence  of  a  magnetic  field,' 
the  electron  drift  velocity  with  respect  to  the  ions  has  only  to  exceed  the 
ion  acoustic  velocity  Cg  to  generate  unstable  waves  traveling  perpendicular 
to  the  magnetic  field.  We  will  consider  the  nonlinear  evolution  of  the 
Farley-Buneman  instability  in  a  low  0,  weakly  ionized,  convecting  plasma 
which  is  subjected  to  a  magnetic  field  Bx,  an  electric  field  Eqz  and  a 
density  gradient  (3nQ/Sz)z.  Differences  in  the  collision  frequencies 
(Vf/Gi  &  1,  ^e/Ge  «  1)  of  the  ions  and  electrons  with  the  background  neu¬ 
tral  gas  results  in  the  formation  of  a  cross  field  current  -JQy  from  the 

Vj  -  E  x  */ B2  electron  drift.  For  weak  currents  Jq  long  wavelength  field 
o 

aligned  fluctuations  in  density  6n  «  exp£i(ky-u)t)]  have  been  found  by 

Simon3  and  Hoh3  to  be  linearly  unstable  when  E  *Vn  >0.  In  the  nonlinear 

~o  o 

regime  this  x  B  gradient  drift  instability  evolves  into  an  isotropic  two- 
dimensional  strongly  turbulent  state  in  the  plane  perpendicular  to  the  mag¬ 
netic  field  as  shown  previously.4  For  stronger  currents  such  that  V  >  C 
the  Farley-Buneman  instability  will  develop  at  shorter  wavelengths.  These 
long  and  short  wavelength  modes  can  coexist  simultaneously  (see  Fig.  1)  with 
the  former  usually  occurring  before  the  latter.  Previous  studies  of  the 
nonlinear  evolution  and  saturation  of  these  short  wavelength  Farley- 
Buneman  modes  have  invoked  quasilinear  effects,5  resonance  broadening,6  and 

mode  coupling.7  These  works  have  neglected  the  effects  of  the  strong  large 
scale  background  E  x  ^  turbulence.  In  this  Letter  we  show  that  the  long 

wavelength  £  x  £  turbulence  can  stabilize  the  short  wavelength  high  fre¬ 
quency  Farley-Buneman  instability.  Although  the  following  discussion  is 
Manuscript  submitted  May  12.  1981. 
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applicable  to  any  weakly  ionized  low  3  current  carrying  plasma  convecting  in 
regions  of  E  x  B  turbulence  it  has  direct  bearing  on  density  irregularities 
in  the  equatorial  electro jet  ionospheric  plasma. 

The  basic  equations  for  the  electron  (Ng)  and  ion  (N^)  fluids  in  a  low 
3,  werkly  ionized,  collisional  plasma  can  be  written  ' 

oN/dt  +  V*NV  =  0  (1) 

eN(E  +  V  x  B)  +  TVN  +  Nm  v  V  =  0  (2) 

~  ~e  ~  e  e~e 


3N/3t  +  V.w  =  0 
~i 


(3) 


(d/dt  +  V .  •7)JJo.V.  =  eN(E  +  V.  x  B)  -  TT>J  -  Nm.v.V. 

~l  i~i  ~  ~  l  i~i 


w 


V*(j.  +  J  )  -  0 

~i  ~e 


(5) 


where  we  have  assumed  quasineutrality  (Ng  ^  =“  >') ,  isothermality  (Te  T^) , 
electrostatic  fluctuations  (£  =  -7$),  and  neglected  electron  inertia. 
Linearizing  equations  ( 1) — ( 5 )  with  N  =  nQ  +  n,  etc.  and  assuming  fluctua¬ 
tions  of  the  form  n,  9,  v^,  vi  <=  exp{i£kvy  +  kzz  -  (u)kr  +  iy(k)t]}with 
k«JB  ■  0  and  k^L  »  1  we  find8  for  the  frequency  and  growth  rate  in  the  ion 
frame 

^kr  "  £’V(1  +  ^  (6) 

/■w 

v<k)  -  C*/(l  +  ♦)]  {  (ne/ve)(Vd/L)cos2e  +  ^r/vt  -  k2cs2/vi}  (7) 

where  \>  m  L  *  *  (l/nQ)(3n  /sz),  cg2  ■  2kT/ra^  and  0  is  the  angle 

defined  by  jk  and  .  From  the  expression  for  the  growth  rate  y(k)  in 


2 


eq.  (7),  we  note  that  at  low  frequencies  (long  wavelengths  and  weak  currents) 
such  that  u-’./v  <  (Cl  /v)  (1/kL)  the  E  x  B  gradient  drift  tern  will  dominate 

K  1  c  6  ^ 

with  all  modes  with  k  <  kc  =  Q  (O^/v^) (V^/L) (v^/cg2) I*5  cos0  unstable.  At 

higher  frequencies  (short  wavelengths  and  strong  currents)  Farley-Buneman 

modes  will  become  unstable  if  >  cg  but  with  no  critical  wavelength.  In  / 

the  equatorial  electrojet  plasma  V,  (t)  is  time  dependent  varying  from 

d 

V  <  c  to  V,  ^  c  over  a  time  interval  At  »  y(k)  1 .  As  a  result  the 
Farley-Buneman  instability  will  be  excited  in  strong  JE  x  J5  turbulence.  Since 
these  waves  are  nondispersive  they  will  interact  strongly.  However  as  their 
amplitude  increases,  they  do  not  steepen  appreciably4  but  are  unstable  to 
perturbations  perpendicular  to  their  propagation.  It  must  be  noted  that  this 
fluid  approximation  is  valid  for  tu  v^.  At  higher  frequencies  kinetic 
effects,  e.g.,  ion  Landau  damping,  will  become  important  and  introduce  a  high 
frequency  cutoff.  For  lower  frequencies  both  fluid  and  kinetic  treatments 
are  identical.9 

By  writing  N(x,t)  =  n  +  L  n(k,uj)exp[,-i(k,x-u)t)  ] ,  <J>(x,t)  =  cp  +£  cp(k,u>) 
~  °kuu~  ~  ~  0  ku>  ~ 

expC~i (k*x-u)t) ] ,  etc.,  expanding  eqs.  ( 1 ) — ( 5 )  in  the  small  parameter 
v^/Ci  <  “i^vi  ^  c  <<:  *  ant^  considering  high  frequencies  uj^  such  that 

”  rw 

uu  /v  >  (Cl  /v  )(l/kL)  10-2  (for  \  3m,  L  6  km,  Q  /v  =  102)  we  find  to 

k  i  e  e  e  e 

second  order  in  n^  (n/nQ  <<:  D 


D(k,w)n  (k,<u)  *  /d2k' du/ V(k,k'  ,u u,u/  )  n(k'  ,uu'  )  n(k-k'  ,u)-u/  ) 
~  ~  J  ~  ~  ~  ~  ~ 


(8) 


where 


_1 


_1  _1 


D(k,uj)  =  uj-k*V  (1  +  p)  -iyw^(l  +  ip)  v  +  iipk2c  2/v 
***  ***  ~d  i  si 

is  the  Farley-Buneman  dielectric  and 
V(k,k'  ,ju .uo7  )  »  -(x  x  k'  *k/k'  2)  [(1  +  y)-1  (v^f.^k'  *Vd  +  k'  *x  x  Vd+ik' 
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By  neglecting  the  nonlinear  term  on  the  right  hand  side  of  eq.  (8)  we 


recover  the  linear  result  from  D(jc,aj^r  +  iyQc))  ■  0  giving  *  jk^V^/d  +  ’!') 

/v  /v 

and  y(k)  *  (iii/1  +  t)(l/vi)  [QS’V^)2  "  k2cs2J' 

Ve  solve  eq.  (8)  for  the  high  frequency  short  wavelength  component  of 

n(k,ou)  by  considering  its  mode  coupling  co  the  low  frequency  long  wavelength 

well  developed  strong  £  x  j5  turbulence.  Let  a  Farley-Buneman  wave  be  denoted 

by  (kj.ujj)  and  a  turbulent  E  x  3  mode  by  (k^.u)^).  Physically,  when  a 

Farley-Buneman  mode  (k  ,ujt)  grows  to  such  a  level  that  it  can  couple  with  the 

^1  X 

E  x  B  turbulence  (kT_,jjT_),  a  beat  wave  component  (kT  db  k _ _ ,  =  jjxt)  will 

~  ~I  III  '-I  I  I  II 

appear  which  in  turn  can  beat  with  (kTT,o)TT)  to  affect  (kT,a)T).  The  evolu- 

~I  I  II  «'■'!  I 

tion  of  n(k_,u)T)  can  then  be  written 
~I  I 

DI(k_,uuT)nI(k  ,uuT)  *  f d2k'  dm  V  (kT  ,k'  ,mT,m  )nII(k/  ,u/  )nX  II(kT-k'  ,uj,-u/  ) 

(9) 

The  beat  wave  (kj-k^  .iu^-id^)  evolves  according  to: 

D1  II(k  -k,uj  -uj)nI  II(k  -k,u)T-o)  =  f d2k!  dm1  V(kT-k,k,,o  -u),uu/  Jn^dk'  ,~m  ) 

nI(kT.dT)  (10) 

I 


where  we  have  used  V(k,k.'  ,u>,u/  )  =  VCk.-jc'  ,u>,-u/ ) .  Substituting  eq.  (10)  into 
eq.  (9)  we  find  the  nonlinear  dispersion  relation  to  lowest  order 


DI(£ru'I)nI(J$i-UJI>  *  0 


(ii) 


T  T 

where  D  (k_,u)T)  *  D  (kT,uaT)  +  6D  (kT,ujT).  The  nonlinear  part  of  eq.  (11) 
~l  I  ~l  l  I 

can  be  written 


r  V(k  ,k'  ,uj  ,u))V(kT-k/  ,k  ,ou  -o'  ,u>' )  TT 

I  /,  „  /,?,/,/  ~I  ~  I  ~I  ~  ~  1  ,11,,./  y  , 

(ir‘"i)  “  •  ^  dUJ  - Ii-ii'”"'/ - 7 - 1  (~  > 

*'  D  (k_-k  ,mT-’x>) 


(11a) 
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r 


b, 

i 
! 


» 


h* 1 


'  i 


where  I^^(k,u))  *  <  n**(k,u>)n**(-jk,-u))/n  2  >  is  the  power  spectrum  of  the 
£  x  _B  turbulence.  It  has  been  previously  shown10  that  III(k,uj)  can  be  calcu¬ 
lated  using  the  direct-interaction  approximation  of  Kraichnan11 


|u>-u»(k)  +  rIX(Jk.ou)  1  2III(k,do)  -  i  J", dV  dw  |  w(k,k7  )|2iri(k7  ,u>7  )IIi:(k-k7  ,iW  ) 

(12) 

TT  r  w(k,k-k/  )w(k-k/  ,k)I1I(k/  ,«/  ) 

T  (k.uj)  -  -  /d2k7du/  - — -  (13) 

^  uu-sju7  -uj  (k-k7  )+r  (k-k7  ,u>-u/  ) 

M  A#  AW  AW 


with  w(k,k/  )  *  V(k,k7  )  +  V(k,k-k7  )  and  uu(k)  ■  u>,  +  iy(k)  where 

-  ~  ~  ~  ~  ~  ~  kr  ~ 

AW 

Y  00  =  (t///l  +  '(/)  f(^g/ve)  (V^/L)cos20  -Jc2cg2/v^].  In  eq.  (12)  r^^(k,u>)  is  the 
self-damping  of  the  long  wavelength  E  x  B  fluctuations  (kTT,iuTT).  We  now  pro- 

AW  AW  /N'II  11 

ceed  to  solve  equations  (11)-(13). 

First,  we  note  that  since  the  interacting  waves  considered  here  are  non- 

I-II 

dispersive  (see  eq.  (6))  the  beat  wave  dielectric  D  (k  -kTT ,uu_-ujtt)  **■  0  and 

~I  ~LL  I  LL 

the  right  hand  side  of  eq.  (11a)  diverges.  As  a  result  we  replace 
D1  **(k  -k7  ,iu  -u/  )  by  its  renormalized  value  (k-k7  ,uj  -u/  )  + 

6D1  ^(k^-k7  .lUj-uj7  )  6D1  (k^-k7  ,ujj-u/  )  in  eq.  (11a).  In  order  to  study  the 

nonlinear  saturation  we  solve  eq.  (11a)  for  Im  6D1  =  -ir(kT  ,uuT(kT) ) .  From 

Awl  I  AVI 

previous  studies4’10  the  steady  state  solution  of  eqs.  (12)— (13)  can  be 
written 


ITI(k,uu)  »  III(k)(21r)'^(rII(k))“1exp[-(u>^(k))2/2(rII(k))2]  (14) 

where  iu(k)  =  Jk»V^/ (1  +  ,  ?ri(k)  *  3.4  n  ls(vi/Qi)k2V£l(III(k))  ,  and 

I^^(k)  *  Ik-n  is  isotropic4’12  with  n  ®*  3-4  while13  I  «  v”,  m  “■  2.  Substitu- 

a 

ting  eq.  (14)  into  (11a)  and  assuming  that  k7  *•  k^  <  k  s  k,  u/  “  < 

uu  s  gu  *^(k„-k7  ,ujt-u/  )  °*  D*(k  ,uuT)  +  0(k7  /kT)  *•  6D*(kT,u)T)  we  find 

1  '"'l  ^  I  *^1  1  1  Awl  1 
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-  «***»  «*‘l 


[r (k)  ] 2  »  J i2k'  v (k , k' ) v (k-k'  , k' )  i11  (k' ) 

“■  (vi/r*i)2^2k/  (x  X  k'  *k) 2  (k'  J-4  III(k/  ) 

=>  (vi/ni)2k2Vd2  JdQ'  k'  dk'  sin2 (Q— 07  )cos20'  Ii:C(k'  )  (15) 


where  r(k)  >>v(k)  has  been  assumed  and  9'  and  9  are  the  angles  made  by  k'  and 
k,  respectively,  with  V..  In  evaluating  the  quantity  V(k,k')  V  (k-k',k) 

^  r*~( 1  ^  <v 

in  eq.  (15)  we  have  kept  only  terms  proportional  to  (v^/fi^)2  where 

>>1.  Experimental  studies14  as  shown  in  Figure  2(a)  indicate  that  T (k) 
is  approximately  independent  of  angle  9.  This  allows  the  replacement  of  T(k) 
in  eq.  (15)  by  its  angle  averaged  result  (2rr)  l^'dQF(k)=:V(k)  giving 


F(k)  *  (vi/Qi)(j)  kVd  <  |n/no|2II  > 


(16) 


where  <  |  n/no|  * JdQ  k  dk  IIX(k)  =  2ir  Jdk  k  III(k).  Using  VjA..  =“  22, 

2i/k  =“  5m,  Vd  =  4  x  102m/sec,  cg  =“  3.6  x  102m/sec,  and  <  fn/nj2^  =“  0.01, 
eq.  (16)  gives  for  the  spectral  width  r(k)/2Ti=“  10  Hz  (cf.  Fig.  1) .  Figure  2(b) 
shows  that  the  scaling  of  T(k)  «  k  from  eq .  (16)  is  in  reasonable  agreement 
with  experimental  results  * 3 ’ 14  which  indicate  that  the  spectral  broadening 
a  F  (k)  «  k2  ‘  7 . 

The  nonlinear  dispersion  relation  in  eq.  (11)  for  the  short 
wavelength  Farlev-Buneman  instability  in  the  long  wavelength  E  x  B  turbulent 
background  can  then  be  written  explicitly 

14  "  +  ^~1_  i^2(l  +  •4/)"1v[1+  i!i/k2cs2/vi  +  ii(vi/f.i)kVd<|n/no|2II>  -  0. 

(17) 

Separating  uu  *  uu^  +  iy(k)  we  find 

y(k)  -  [*/(l  +  il>)]  |(k*Vd)2/vi  -  k2cg2/v1f  -  j(vi/ai)kVd<Sn/no|2II>  (17a) 
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r- 


For  complete  stabilization  of  the  fastest  growing  linear  mode  (8  0)  the  long 

wavelength  background  fluctuation  level 

<  jn/n  j2  £  2  v  0  (1  +  ip)  1  (k V./v.  -  k2c  2/v  kV  ). 

O  II  66  Q  1  s  la 

with  kV,/v,  0.1,  k2c  2/v.kV,  “"0.05  we  find  n/n  ^  0.002  which  is  consis- 
d  i  sid  o 

tent  with  available  experimental  estimates  . 
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Doppler  Shift  (Hz) 

Fig.  1  —  Simultaneous  equatorial  electrojet  ionospheric  plasma  density  fluctuation  power 
spectra  (vertical  axis  I(k  ,o>)  vs.  frequency  u  in  Hz  (horizontal  axis)  from  Ref.  14  for  several 
radar  backscatter  observation  frequencies;  the  radar  frequencies  29.0  MHz,  9.0  MHz, 
6.6  MHz,  5.6  MHz  correspond  to  wavelengths  5.1m,  10.5m,  16.6m,  22.7m,  and  26.7m, 
respectively. 
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Fig.  2  —  (a)  Power  spectrum  width  vs.  wave  vector  angle  made  by  k  and  Vd,  the  electrojet 
electron  drift  velocity,  at  several  wavelengths  from  Ref.  14;  note  approximate  isotropy  in 
angle,  (b)  Power  spectrum  width  vs.  wave  number  k  at  9  -  15°  from  Ref.  14;  solid  line  is 
best  fit  a  k0-'. 
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